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We show that the recent reinterpretation of oxygen isotope efTects in cuprate superconductors by 
D. R. Harshman et al. is mathematically and physically incorrect violating the Anderson theorem 
and the Coulomb law. 
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The doping dependent oxygen isotope effect (OIE), a, 
on the critical superconducting temperature Tc (for re- 
cent reviews see Ref.^) and the substantial OIE on the 
carrier mass^, Um* , provide direct evidence for a signifi- 
cant electron-phonon interaction (EPI) in cuprate super- 
conductors. High resolution angle- resolved photoeniis- 
sion spectroscopy (ARPES)^ provides further evidence 
for the strong EPI^ apparently with c-axis-polarised 
optical phonons. These results along with optical^'^'^ 
and neutron scatteringS>S spectroscopies unambiguously 
show that lattice vibrations play a significant but un- 
conventional role in high-temperature superconductivity. 
The interpretation of the optical spectra of high-Tc ma- 
terials as due to many-polaron absorption^" strength- 
ens the viewii that the Frohlich EPI is important in 
those structures. Operating together with a short-range 
deformation potential and molecular- type (e.g., Jahn- 
Teller— ) EPIs, the Frohlich EPI can readily overcome 
the Coulomb repulsion at a short distance of about the 
lattice constant for electrons to form real-space intersite 
bipolarons or Cooper pairs depending on doping^"^. 

Despite all these remarkable and well-done experi- 
ments that lead to the consistent conclusion about the 
important role of EPI in high-temperature superconduc- 
tors, Harshman et al.^'^ have recently claimed that the 
observed large OIE is caused by a disorder-induced pair- 
breaking rather than by strong electron-phonon coupling 
and/or polaronic effects. Based on their reinterpretation 
of OIE, they conclude that EPI is allegedly too weak 
to explain high Tc in all the high-Tc materials. Here 
we show that the reinterpretation of OIE^^ is internally 
inconsistent being at odds with a couple of fundamen- 
tal physical laws. More specifically we show that the 
reinterpretation stems from a mathematically incorrect 
formulism. 

Given the added claim that the pairing symmetry is 
nodeless s- wave, the authors in Reft^ have assumed that 
the variation of T^ with doping is determined by the "uni- 
versal" relation, 

ln(Tco/rc) - a[ijj{\/2 + 1/Tct) - VXl/2)], (1) 

which was originally derived by Abrikosov and Gor'kov™ 
with the coefficient a = 1 to describe the pair-breaking ef- 
fect by magnetic impurities in conventional s-wave BCS 
superconductors. Here Tco is the critical temperature 



of optimally doped compounds in the absence of pair- 
breaking, i){x) is the digamma function, and r = AnTtr 
is proportional to the transport relaxation time, rt^, due 
to impurities, which are thought to be responsible for 
the suppression of TcO (we take h = ks = ^ here and fur- 
ther). Since nonmagnetic disorder in cuprate supercon- 
ductors also often depresses Tco, Harshman et alM- have 
erroneously relaxed the requirement of magnetic impuri- 
ties applying Eq. ([1]) to nonmagnetic impurities with the 
same coefficient a = 1. 

In fact, the coefficient a in Eq.([T]) strongly depends on 
the pairing symmetrjslfi. as analyzed in detail by Fehren- 
bacher and Norman^S. For nonmagnetic impurities a — 1 
holds only for a d-wave (DW) or g-wave (GW) super- 
conductor with a zero average gap, while this coefficient 
is significantly smaller in an anisotropic s-wave (ASW) 
superconductor^". When the BCS gap is isotropic, the 
familiar "Anderson theorem", Tc = Tco, is satisfied^ii^S 
because a — 0. But even in the extreme case of a highly 
anisotropic ASW superconductor with the same nodal 
structure as in the DW superconductor the effect of non- 
magnetic impurities on their properties remains qualita- 
tively different, although the two states are indistinguish- 
able in phase-insensitive experiments^. In particular we 
show here that the pair-breaking OIE enhancement is 
negligibly small based on any s-wave gap function that 
does not change sign with angle, contrary to Ref.^''. 

On the other hand, the effect of magnetic impurities 
in an ASW superconductor, or the effect of nonmagnetic 
impurities (or disorder) in a DW or GW superconduc- 
tor can cause a significant enhancement of the isotope 
effects on both Tc and the penetration deptbi^ii^. Two 
different groups"'^^'^^ have consistently shown that the iso- 
tope effects on Tc and the penetration depth are almost 
proportional to each other provided that the strong pair- 
breaking effect exists. These theoretical models may be 
able to explain the observed large oxygen-isotope effects 
on both Tc and the penetration depth in underdoped 
cuprates if the scattering rate were large enough. How- 
ever, these models cannot consistently explain the negli- 
gibly small OIE on Tc but a large OIE on the penetration 
depth in optimally doped cuprates^. 

Differentiating Eq.((T]) with respect to the ion mass, 
M, one can express OIE, a — —dlnTc/dlnAI, in 
terms of the OIE observed in optimal compounds, ao = 



a / ao 
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FIG. 1: Pair-breaking enhancement of the oxygen isotope ef- 
fect, a/ao, in the d-wave (DW) superconductor (left panel) 
and in an anisotropic s-wave (ASW) superconductor with the 
same nodal structure (right panel) as a function of the pair- 
breaking parameter, 1/r. Symbols (right panel) represent the 
experimental data for Zn-doped YBa2Cu307_a used in Ref.~. 



-dlnTco/dlnM, 
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l-aii)i{l/2 + l/T,T)/Ty 



(2) 



where '4'i{x) = dip{x)/dx is the trigamma function, if r is 
independent of M. As shown in Fig. 1 (right panel) using 
Eq.(l2]) the maximum OIE enhancement is about 30% or 
less even in the extreme case of the ASW superconductor 
with the same nodal structure as in the DW superconduc- 
tor, where the enhancement is huge, about several hun- 
dred percent or more. Fig. 1 (left panel). For a nodeless 
s-wave gap, hypothesized in Ref.^*^, there is practically no 
enhancement at all. Mathematically the difference comes 
from the different numerical coefficients in Eq.([T]): a — 1 
for DW, a = 1/4 for the extreme ASW^S, and a < 1/4 for 
a nodeless gap. Physically the difference comes from the 
non- vanishing, impurity-induced, off-diagonal self energy 
in the ASW state, which is absent in the DW stat o^^'^° . 
As a result the "pair-breaking" reinterpretation of OIE 
by Harshman et al^ with the nodeless pairing symme- 
try turns out to be incompatible with the experimental 
data. The experimental OIE, Fig. 1 (left panel), is more 
than one order of magnitude larger than the predicted 
OIE when the correct equation is applied. Fig. 1 (right 
panel). 

One can also rule out the pair-breaking explanation 
of OIBi^ even in the case of the DW order parameter, 
in particular in Pr substituted YBa2Cu3 07_j^ (YBCO), 
although there is apparently good agreement with the 
experiment in the case of Zn-doped YBCO, as seen in 
Fig. 1 (left panel). Since Zn doping induces a magnetic 
moment of about 0.8 /is per Zn, the data might be con- 
sistent with the magnetic pair-breaking effect in the case 
of an s-wave symmetry. But for YBCO with oxygen 
vacancies or substituted by trivalent elements for Ba, 
no magnetic moments and disorder are induced in the 
Cu02 planes so that the impurity scattering rate may 
increase only slightly. In fact the low-temperature co- 
herence length in cuprate superconductors is very small, 



5o = Q-18vf/Tc < 2nm, while the mean free path, I = 
vpTtr, is about 10 nm or larger as follows from resistivity 
and recent quantum magneto-oscillation measurements 
in the underdoped YBa2Cu306.5^^ {vp is the Fermi ve- 
locity). Using these data one obtains 1/TcT < 0.1, which 
is too small to account for the observed enhancement 
of OIE with any gap symmetry as seen from Fig. 1, or 
for the doping dependence of the magnetic penetration 
depth, contrary to Ref.^**. 

If however, in spite of the above estimate, Pr substi- 
tution might lead to a pair-breaking parameter, 1/TcT > 
0.5, large enough to explain the enhancement of OIE on 
Tc, one should expect a similar enhancement of OIE on 
the penetration depth because the magnitudes of the en- 
hancement in the isotope effects on Tc and the pene- 
tration depth are nearly proportional to each other— li^. 
Nevertheless, OIE on the penetration depth is nearly con- 
stant from the optimally doped sample to the substituted 
samples with a large amount of Pr^**, which is inconsis- 
tent with the theoretical predictio n^^'^^ . Claiming the 
opposite, Harshman et al. have made further mistakes in 
their derivation of the penetration depth, \ab (Eq. (8) in 
RefJ^). They have applied the conventional correction 
factor l-l-^o/^ due to impurity scattering, which actually 
yields 



\l,{Tc)^\l,{T,,)[l + QMa/T,_ 



(3) 



where a = l/2Tt,.. Eq.dS]) differs from Eq. (8) in RefJ^ 
with Tco instead of T^ in the second term inside the 
square brackets. Clearly using Eq.(I3]) instead of the in- 
correct Eq. (8) in Ref.^'* one obtains an enhancement of 
OIE on \ab similar to that on T,. contrary to the erro- 
neous claim of Refpi^. Moreover the penetration-depth 
formula Eq.([3]) is valid only for nonmagnetic impurities 
in s-wave superconductors that do not suppress T^. If one 
assumes that nonmagnetic impurities can suppress T^, 
one should consistently use a formula for the penetration 
depth, which is also associated with the pair-breaking 
effect. 

We would like to emphasize here that, since the pair- 
breaking effect in optimally doped cuprates is negligibly 
small and the carrier concentrations of the two oxygen- 
isotope samples have been consistently proved to be 
the same within ±0.0002 per Cu^i^i, the observed large 
oxygen-isotope effect on the penetration depth must be 
caused by the large oxygen-isotope effect on the supercar- 
rier mass. The origin of this unconventional isotope effect 
should arise from strong EPI that causes the breakdown 
of the Migdal approximation. Indeed a model based on 
(bi)polarons^'^ accounts naturally for both OIEs, a and 
Urn* ■ There is a qualitative difference between ordinary 
metallic and polaronic conductors. The renormalized ef- 
fective mass of electrons is independent of the ion mass 
M in ordinary metals (where the Migdal adiabatic ap- 
proximation is believed to be valid), because the EPI 
constant A = Ep/D does not depend on the isotope mass 
{D is the electron bandwidth in a rigid lattice). How- 
ever, when electrons form polarons dressed by lattice dis- 



tortions, their effective mass m* depends on M through 
m* = m exp(7i?p/w), where m is the band mass in a rigid 
lattice and 7 < 1 is a numerical coefficient depending on 
the EPI range. Here the phonon frequency, w, depends 
on the ion mass, so that there is a large polaronic iso- 
tope effect on the carrier mass with the carrier mass iso- 
tope exponent am* = dh\m* /d\nM = (1/2) ln(7n*/TO) 
as observed^, in contrast to the zero isotope effect in or- 
dinary metals. Importantly a™* is related to the critical 
temperature isotope exponent, a, of a (bi)polaronic su- 
perconductor as Of = Um* [1 — ("z/m*)/(A — //c)], where ji^ 
is the Coulomb pseudo-potentiaP^. Contrary to another 
misleading claim by Harshman et a/J^, the latter expres- 
sion accounts for different doping dependencies of a and 
am* as well as for a small value of a compared with «.,„. 
in optimally doped samples, where the electron-phonon 
coupling constant A approaches from above the Coulomb 
pseudopotential jjl^. Similarly, the unconventional iso- 
tope effects^ were also explained by polaron formation 
stemming from the coupling to the particular quadrupo- 
lar (5(2)-type phonon mode in the framework of a multi- 
band polaron model-^. 



Finally the claim by Harshman et ali^ that EPI is 
weak in high-Tc superconductors compared with the 
Coulomb coupling between carriers in buffer and Cu02 
layers contradicts the Coulomb law. EPI with c— axis 
polarized optical phonons is virtually unscreened since 
the upper limit for the out-of-plane plasmon frequency^^ 
(< 200 cm~^) is well below the characteristic frequency 
of optical phonons, ui « 400 -^ 1000 cm~^ in all cuprate 
superconductors. As the result of poor screening the 
magnitude of an effective attraction between carriers, 
Vphir) = — e^(e^^ — ep ^)/r, induced by the Frohlich 
EPI, is essentially the same as the Coulomb repulsion, 
Vc{r) = e^/eooT, both of the order of 1 eV, as directly 
confirmed by a huge difference in the static, eo 3> 1, and 
high-frequency, £00 « 4 -^ 5, dielectric constants of these 
ionic crystals?^. 

To summarize we have shown that the conclusions 
by Harshman et alM- are mathematically erroneous and 
physically at odds with the fundamental Anderson theo- 
rem and the Coulomb law. We thank Annette Bussmann- 
Holder for calling our attention to Refi^^ and illuminating 
discussions. 
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